A Projection Method Based on Extended Krylov Subspaces for Solving Sylvester Equations
نویسندگان
چکیده
In this paper we study numerical methods for solving Sylvester matrix equations of the form AX +XB +CD = 0. A new projection method is proposed. The union of Krylov subspaces in A and its inverse and the union of Krylov subspaces in B and its inverse are used as the right and left projection subspaces, respectively. The Arnoldi-like process for constructing the orthonormal basis of the projection subspaces is outlined. We show that the approximate solution is an exact solution of a perturbed Sylvester matrix equation. Moreover, exact expression for the norm of residual is derived and results on finite termination and convergence are presented. Some numerical examples are presented to illustrate the effectiveness of the proposed method. Keywords—Arnoldi process; Krylov subspace; Iterative method; Sylvester equation; Dissipative matrix
منابع مشابه
On the numerical solution of generalized Sylvester matrix equations
The global FOM and GMRES algorithms are among the effective methods to solve Sylvester matrix equations. In this paper, we study these algorithms in the case that the coefficient matrices are real symmetric (real symmetric positive definite) and extract two CG-type algorithms for solving generalized Sylvester matrix equations. The proposed methods are iterative projection metho...
متن کاملOn the Numerical Solution of Generalized Sylvester Matrix Equations
The global FOM and GMRES algorithms are among the effective methods to solve Sylvester matrix equations. In this paper, we study these algorithms in the case that the coefficient matrices are real symmetric (real symmetric positive definite) and extract two CG-type algorithms for solving generalized Sylvester matrix equations. The proposed methods are iterative projection methods onto matrix Kr...
متن کاملTheoretical results on the global GMRES method for solving generalized Sylvester matrix equations
The global generalized minimum residual (Gl-GMRES) method is examined for solving the generalized Sylvester matrix equation [sumlimits_{i = 1}^q {A_i } XB_i = C.] Some new theoretical results are elaborated for the proposed method by employing the Schur complement. These results can be exploited to establish new convergence properties of the Gl-GMRES method for solving genera...
متن کاملSylvester Equations and the Factorization of the Error System in Krylov Subspace Methods
This paper presents a factorization of the error system that arises in model reduction of linear time invariant systems by Krylov subspace methods. The factorization is introduced for reduced models that match moments and/or Markov parameters of the original system with multiple inputs and outputs. Furthermore, dual results are given for the reduction with input and output Krylov subspaces. To ...
متن کاملPreconditioned Galerkin and minimal residual methods for solving Sylvester equations
This paper presents preconditioned Galerkin and minimal residual algorithms for the solution of Sylvester equations AX XB = C. Given two good preconditioner matricesM and N for matrices A and B, respectively, we solve the Sylvester equations MAXN MXBN =MCN. The algorithms use the Arnoldi process to generate orthonormal bases of certain Krylov subspaces and simultaneously reduce the order of Syl...
متن کامل